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Abstract 

We calculate the two-point correlation function and magnetic susceptibility in 
the anisotropic 2D Ising model on a lattice with one infinite and the other finite 
dimension, along which periodic boundary conditions are imposed. Using exact ex- 
pressions for a part of lattice form factors, we propose the formulas for arbitrary 
spin matrix elements, thus providing a possibility to compute all multipoint corre- 
lation functions in the anisotropic Ising model on cylindrical and toroidal lattices. 
The scaling limit of the corresponding expressions is also analyzed. 

1 The model 

Two-point correlation function and magnetic susceptibility of the isotropic Ising model on 
a cylinder were calculated in [Tl|2]. Analogous results can be obtained in the anisotropic 
case as well. The complications due to the presence of two different coupling constants can 
be overcome, since the matrix whose determinant yields the correlation function continues 
to have a Toeplitz form. The computation idea has already been presented in |pLj, and 
below we will often use the results of that paper, omitting detailed calculations. 

The calculation of multipoint correlation functions can be reduced to the problem of 
finding Ising spin matrix elements (the so-called form factors) in the orthonormal basis 
of transfer matrix eigenstates. Although neither the first nor the second problem has 
been solved on a finite lattice, the language of matrix elements turns out to be more 
convenient when constructing the corresponding expressions. A part of form factors can 
be found from the formulas for the two-point correlation function. Since the structure of 
these expressions is relatively simple, we generalize them to the case of arbitrary matrix 
elements. 

The hamiltonian of the anisotropic Ising model on a rectangular lattice is defined as 

H[cr] = - 5^ a(r) ( J,.V.. + JyVy) a{v), 

T 

where the two-dimensional vector r = (r^,ry) labels the lattice sites: r^. = 1, 2, . . . , M, 
Vy = 1, 2, . . . , A^; the spins cr(r) take on the values ±1. The parameters and Jy 



1 



determine the coupling energies of adjacent spins in the horizontal and vertical direction. 
The operators of shifts by one lattice site, and Vy, are given by 

'^xO-{r^,ry) = a{r^ + l,ry), Vy(j{r^,ry) = cr{r^,ry + 1), 

where for periodic boundary conditions one has 

(vf ')"' = (vf ))" = 1, 

and for antiperiodic ones 

(Vf^))^'= (Vf^))^ = -1. 

If the lattice is periodic in both directions, the partition function of the model at the 
temperature 



can be written as a sum of four terms 

each of them being proportional to the pfaffian of the operator 

/ 1 1 + 1 \ 

D = 



4 1 1 + 



-l-t,V_, -1 1 

V -1 -l-tyV-y -1 / 

Q = (2 cosh cosh Kyf^^ ■ Pf Z) , 

with different boundary conditions for V^., V^,. Here, we have introduced the dimension- 
less parameters 

= jSJx, Ky = pjy, tx = tanhi^^j., ty = tanh.K'j^. 

One can verify that when the torus degenerates into a cylinder (M ^ A^), the partition 
function is determined by the antiperiodic term: Z = Q^^^'^^\ 

Let us consider the two-point correlation function in the Ising model on the cylinder 

(a(0, 0)a(r., r,)) = Z'^ J] a(0, 0)a(r., r,)e-^^W = = (2) 

where Zdef denotes the partition function of the Ising model with a defect: the coupling 
parameters Kx, Ky should be replaced by Kx — iTT/2, Ky—m/2 along a path that connects 
the correlating spins (see Fig. 1 in [1] , where the numbering of lattice sites and the locations 
of correlating spins were described; the bold line in the figure is the path along which the 
couplings are modified). The exponent P in the formula ([2]) is equal to the number of steps 
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along the defect line: P = + Vy in the case of a shortest path. When the correlating 
spins are located along a line parallel to the cylinder axis (i. e. = 0), the ratio of 
pfaffians in the right hand side of ([2]) can be expressed in terms of the determinant of a 
Toeplitz matrix 

(a(0,0)a(r,,0)) = det/l, (3) 

~ J 2TrN ^ (1 + tl) (1 + tl) - 2t, (1 - tl) cosp - 2ty (1 - tl) cosg' 

(4) 

whose size \rx\ x \rx\ is determined by the distance between the correlating spins. Here 
and below the superscripts (NS) and (R) in sums and products imply that the cor- 
responding operations are performed with respect to Neveu-Schwartz {q = ^ (j + |) , 
j = 0, 1, . . . , — 1) or Ramond (g = ^ j, j = 0, 1, . . . , — 1) values of quasimomenta. 
Our goal is to transform into a representation with an explicit dependence on the 

distance. 



2 Ferromagnetic phase 

In the translationally invariant case the pfaffians Pf D can be easily calculated, 

Q = J^(cosh 2Kx cosh 2Ky — sinh 2Kx cos q — sinh 2Ky cosp)^/^. (5) 

The product over any of the two quasimomentum components in ([5]) can be found in 
an explicit form. For instance, the term in ([1]), which corresponds to periodic boundary 
conditions along the x axis and antiperiodic ones along the y axis, may be written as 

QiR,NS) ^ (2sinh2i^,)^^/2-Q(«)g7V,(,)/2^^^g-7V,(,))^ 

= {2smh2K,)^^/^l[^'"\^^^Py\l-e-^^^P^), (6) 

q 

where the functions 7(g) and 7(p) are determined by the relations 

(l+tg.)(l+tg) ty{l-tl) 



'^'^^(p^= 2tyii-ti) — ^r^''''^ 



and the conditions 7(g), 7(p) > 0. 



We will call the domain of values of K^, Ky, where 

sinh 2/^3. sinh2Ky > 1 
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(and hence 7(0) = lntx + 2Ky) the ferromagnetic region of parameters. Notice that in this 
case the numerator of the integrand in (jlj) can be represented in the following factorized 
form in terms of the function 7(g): 

2tx (1 + tl) + {tl - 1) {e-'P + tyP) = (1 - tl)e"'^°\l - e-^We'P)(l - e-^^^^e-'^^'). 

Then, using the identity 



N ^ 



(NS) 



tanh^ 



cosh 9 — cos q sinh 9 



(9) 



one may compute the sum over the discrete Neveu-Schwartz spectrum in (jl]). As a result, 
we can express matrix elements Akk' in terms of contour integrals 



Akk' 



2m 



-z''-'A{z), 
z 



A{z) 



where 



(l_e-7W^)(l_e-T(o)z-i) 
_(l_e-7(o)2)(l_e-7W2-i)_ 

iV7(p) 



T{z) 



T{z) = tanh ■ 



z = e'P. 



(10) 



(12) 



Recall that in order to calculate the determinant |A| by the Wiener-Hopf method [3] 
one has to represent the kernel A{z) in a factorized form 



A{z) = P{z)Q{z-'), 
where the functions P{z) and Q{z) are analytic inside the unit circle. Setting 

P(z) 



W{z) 



the determinant can be written as follows: 

00 

where 



1=0 



h{r,) = r, [lnP(0) + lnQ(0)] + 



2ni 



d 

dz\nQ{z'^)—\nP{z) 



(13) 



(14) 



(15) 



(16) 



1^1=1 



m\{27riy 



21 i-i I I 

Uidz,zl^)Yl n [iz2^-l-Z2j-ifiz2.-Z2,y] YlWiz2^-l] 
1=1 i=lj=i+l j=x 



n 11(1-^2,-1^2,)' 

i=l j=l 



I 

Y 

i=l 



n w{z,^) 



(17) 
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We will see that the sum (llSp contains a finite number of terms, i. e. there exists Iq such 
that g2i = for all I > Iq. 

Let us now calculate the integral (|T71) . To do this, we first rewrite the function T{z) 
defined by (fT2l) in the factorized form: 



T{z) = tanh 



Y[^^\cosh^{p) — cosg) 
<? 

n'-^'^^(cosh7(p) — cosg) 



L <? 



q 

|-[(A^s)g7(<?)(i - e-^('?)2)(l - e-^W^-i) 

L 9 

Here, we have introduced the notation 



Next let us define 



n(-)(i-e 



-7(9) . 



]^{^^)(l_e-7(</)^) 



Qriz) 



1 - e-^Wz' 
1 - e-^(o)^ 



and write as follows: 



W{z) 



Qo{z) 



1 - e-^^^h 
1 _ e-TW; 



4.(1 -t^) coth^^SM 
^ (cosh 7(7r) - cos p) ^ 



7(0)-7(7r) 



-2A- 



tyil-tl) 



where z = e*^. Taking into account that 

1 ty{l - t 



coth™ 



sinh7(p) Nt,{l-tl)^ 

y q 



one can write W{z) as 

W{z) = (cosh 7(71) — C0SJ9) Qt{z) e 



(R) 



sinh7(p) 



1 



cosh 7(g) — cosp 



-y 

AT ^ 



(R) 



(1^ 



(19) 



N ^ cosh 7(g) — cosp 



Analogous computation forW ^{z ^) gives 



(cosh 7(0) — cosp) Qt{z) e 



y(^)-7(0) 1 -sr^ifi) 



-y 



N ^ cosh 7(g) — cosp 



Using the last two formulas, we can calculate the integral (jlTl) . Let us interchange the 
order of integration and summation over discrete quasimomenta. Then all the singularities 
inside the integration contours are exhausted by a finite number of poles at the points 

27]-/- 



corresponding to the Ramond values of quasimomentum. The integral ( fT7|) can then be 
computed by residues, and the result reads 



g-2«/A (R) ^ g-r:r7{g»)-'?{g») 
^2'(-^)=(^L 11 sinh7(g.) ^^'f'^' 



where 



_ J. 



-7(g)-7('?i)'i 



The function Q2i\(l\ is given by 



QM = c^'u n 



i=l j=i-\-l 

I 

n 

1=1 



-7fei-i) _ ^-7(923-1)^^ ^p-7(i?2i) _ p-l{q2j) 



X 



cosh7(7r) - cosh7(g2j-i) I ( cosh 7(^22) - cosh 7(0) 



21 I I 

Y[ e'^fe) n n (1 - e-'^te»-i)-7(g2j)y 

i=l 1=1 j=l 

l-l I 



(20) 



CI 

22P 



21 r 



2\-\2l 



n n [(cosg2i-i - cosg2i-i) (cosg2i - cosg2j)' 



i=l j=i+l 



21 

Yl sinh 

i<j 



X 



2 7(gi)+7(gj) 



J]^ [(1 + cosg2i-i)(l - cosg2i)] • 



It can be symmetrized with respect to the permutations q2i ^ Q'2j-i- It turns out that 
the result of such a symmetrization coincides with the even part of the function 



ty (1 - tj) 



2/2 21 . 2 qi-qi 
sm 



n 



„;„|,2 7ta)+7fe) ' 
i<j 2 



so that g2i{fx) niay be written as 



92i{rx) 



-2//A 



{2i)\m 



21^ 



{R)e 



-\r:c\-i(qi)-'q(q^) 



smh7(gi) 



^70 = 1, (21) 



Fn[q] = n 



sm 



■x. 7fe)+7(gj) 



(22) 



The formula fl22l) imphes that for 21 > N F2i[q] = 0. By virtue of the relation f|T5l) . we 
then obtain the following representation for the correlation function in the ferromagnetic 
region: 



where 



[N/2] 



(o-(0,0)a(r.,0))(-) = e^^^l/^ J] g^nir.) 



n=0 



^ = |l - (sinh2ir^sinh2i^'y)-2| 

n^'^^n^^'^^sinh^^^^ 

q p 



Y{^^^Y{^^^ sinh n^^'^'^n^'"''^ sinh 

q p q p 



(23) 
(24) 
(25) 



3 Paramagnetic phase 

In the paramagnetic region of parameters, determined by the condition sinh 2Kx sinh 2Ky < 
1, we have the inequality Int^. + 2Ky < 0, and hence 

7(0) = -lnt,-2Ky. 

Because of this, the kernel of the Toeplitz matrix whose determinant yields the correlation 
function should be rearranged: 



A{z) 



(a(0,0)cr(r„0))(+) = \A\, 
(l_e-7W^)(l_e-^W^) ■ 



T{z). 



(26) 



(27) 



_(l_e-7(o);2-i)(l_e-7W^-i)_ 

To calculate the determinant of such a matrix, we must modify the method described in 
the previous section. For the factorized kernel 



A{z) = Piz)Q{z-'), 
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where the functions In P{z) and \nQ{z) are analytic inside the unit circle \z\ < 1, we have 
the following general formula [1] 

1-41= ^ ^^,(0) J^h^iir^), (28) 



1=0 

I 



f^y J. 21+1 n ( Z2iW{z2i 



1^1 *=1 



l-l I I l+l 

n n {z2i-z2jfY[ n (2:2^-1-^2^-1) 

1=1 j=i+l 1=1 J=j+1 

n n (1 - z2.^2,+i)2 



(29) 



i=l j=0 

Here, 1^(2;) and h{rx) are defined by the formulas ( fT4|) and ( fT6l) . Using the transformations 
analogous to those described above, one finds 



zW{z) = 2ze 2 -^/^Ql(^) (cosh 7(0) - cosp)(cosh 7(7r) - cosp) X 

1 v^(fi) 1 



X 



N ^ cosh 7(g) — cosp 



cosh 7(g) — cosp 

Using the last two relations, the integral (12^ can be computed by residues: 



2i+l 



where 



_ 2^^+2/ ^_2Z/A+ 2(0)+2M ^(^^ 2|W ^j^j 



En 



(2/ + i)!iV2m 11 3inh7(g,) ^^^^.^^2 



•^[^] = 22^Wi) n n (cosg2i - cosg2j)^n n (cosg2i-i - cosg2i-i)^ JJ ^^^^ • 

i=\ j=i+\ i=l j=j+l i=l 

We again note that after symmetrization over the permutations q2i g2j+i, the function 
T\q\ coincides with the even part of the function 



2«+l 

n. 2li~ 
sm 

i<j 
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The correlation function in the paramagnetic region of parameters is then given by 



[{N-l)/2] 



-(2Z+1)/A 



92i+i[rx) 



(2/ + l)!Ar2m 



n=0 



Fn[q] = n - 



sm 



sinh 



7fe)+7(gj) 



sinh7(gj 
, F, = l. 



(30) 



i^lmi?], (31) 



(32) 



«<J """" 2 

Comparing the formulas (!2T|) - (!25|) and (!30l) - (l32l) for the correlation function with the 
isotropic case, we see that all the difference consists in the redefinition of the function 
7(g) given by ([7]) and in the appearance of the multipliers 

tyil-tl) 



in the appropriate powers in the function gn{fx)- In the thermodynamic limit N ^ 00, 
these expressions reduce to the known formulas 



4 Magnetic susceptibility 

Since we have two independent coupling parameters in the anisotropic Ising model on the 
cylinder, we can consider two limit cases: 

^0, KyJ^O 

and 

^0, Ky-^ 0. 

The two-dimensional Ising lattice splits into mutually non-interacting closed one-row Ising 
chains of the length in the first case, and into chains of infinite length in the second 
case. Such limit cases provide primitive tests for the form factor representations obtained 
above. To perform such checks, it would also be desirable to have an expression for the 
correlation function {a{0)a{r)) of two spins with arbitrary location on the lattice, i. e. 
for 7^ and Vy 7^ 0. The same problem arises when one tries to write the lattice 
analog of the Lehmann representation for the two-point Green function or to compute 
the magnetic susceptibility: in both cases we have to perform the Fourier transformation 

e*'i'"(o"(0)o"(r)). Meanwhile, the matrix whose determinant gives the correlation func- 
tion has a non-Toeplitz form for Vy 7^ 0, and the method used above to find form factor 
expansions becomes inapplicable. The answer can nevertheless be obtained in this case 
as well if we formulate the problem in terms of spin matrix elements. 

Recall that the 2^-dimensional space in which the Ising transfer matrix T acts can be 
splitted into two invariant subspaces — the Neveu-Schwartz sector (NS) and the Ramond 
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sector (R). The eigenvalues corresponding to the eigenvectors from each subspace are 
given by 



A(-^) = (2sinhirX/^expi{±7(^)±7(|)±...±7(2^-^)}, (33) 

A(«) = (2sinhi^^.)^/^expi|±7(0)±7(|)±...±7(2vr-|)}, (34) 

and we have different selection rules in different regions of parameters of the model: the 
number of minus signs in the NS-sector (l33l) is even in both ferromagnetic and para- 
magnetic region, while the number of minus signs in the R-sector flM|) is even in the 
ferromagnetic region and odd in the paramagnetic one. Note that the expressions fl5^ . 
fl5^ for the tranfer matrix eigenvalues can be obtained from the representation ([T]) for 
the partition function, if we expand the products ([5]): 

i 

ZiR) = l(QiNS,R) _ QiR,R)-^ ^ 5^(Af )^^ 

2 

i 

Because of translation invariance, the transfer matrix can be diagonalized simultane- 
ously with the translation operator. This additional requirement determines the orthonor- 
mal basis of eigenvectors up to permutations. It is convenient to interpret the elements 
of this basis in terms of NS- and R-multiparticle states with appropriate values of quasi- 
momenta: they are generated from the vacua |0)ivs' and \^)r by the action of fermionic 
creation operators. 

Let us consider for definiteness the ferromagnetic region and show how to express the 
two-point correlation function in terms of the matrix elements 

Nsiqi, ■ ■ • ,g2x|'3-(0,0)|pi, . . .,P2l)r 
of Ising spin (the elements NS-NS and R-R vanish due to Z2-symmetry of the model): 

(a(0,0)a(r„r,)) = ^s(0|^(O, 0)a(r,, r,)|0);v5 = 

= ^s(0|a(O,O)r^^a(O,r,)r-^^|0)^s = 

-^^M 

= E 7^ E \NsimoMpi---P2K) 



R\ 



X 



K=0 ^ ' P1---P2K 

, 2K 2K ^ 

X exp ( -r^ -f{pj) + iVy ^Pj) 
V ,=1 ,=1 / 



(35) 



10 



Comparing this expression with the formulas (pT]) - fl23l) . we find spin matrix elements 
between the NS-vacuum and an arbitrary R-eigenstate: 



\NS\ 



i<T(o,o)bi...pz.)ij|' = eeT 



n 



\Nsmh-f{pj) f-^-^sinh 



n 



sm 



2 Pi-Pj 



,2 7(Pi)+7{Pj) 



(36) 

The formula ( l35l) then implies that form factor expansion of the correlation function 
((j(0, 0)a{rx, Ty)) is obtained from (|2Tl) -(|23l) by the substitution 

^-\'rxb{qi)-r){qi) Q-\rx\l{qi)+iryqi-r){qi) 



sinh7(gj 



sinh7(gj 



in the function Qn- An analogous result may be found in the paramagnetic C3iSG clS well. 
Thus we have 



[N/2] 



(a(0,0)a(r,,r,)) 



(-) 



1=0 

m-i)/2] 



(a(0,0)a(r.,r,))W = ^Te-I^^l/^ ^72m(r.,r,), 

1=0 



(37) 
(38) 



where 



-n/A 



tx(i-€). 



q i=i 



-k^l7fe)+«'"Hgj-'?{gj) 
sinh7(gj) 



Fnb] ■ (39) 



As an illustration, consider the "paramagnetic" expansion for Ky = 0. Notice that in 
this limit 

cosh 7(g) = coth2ii'3., 7(g) = = const, 

and, therefore 

^ = ^T = 1, A-i= 77(g) = 0, 



(a(0,0)a(r.,r,)) = 1 5^^^t.)l-le-^'^ = (t.)l-l S 



Ory 



As expected, the spins from different lattice rows are uncorrelated. 

Magnetic susceptibility of the two-dimensional Ising model on a finite M x N lattice 
in zero field may be written as a sum of correlation functions: 



AI-l N-1 



In order to find the susceptibility on the cylinder, we can use the expression 

00 N-l 
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with some precautions f2\. Because of simple structure of the form factor expansions 
fl37|) - fl39l) this sum is easily computed, and we obtain the following expressions for the 
susceptibility in the ferromagnetic and paramagnetic region of parameters: 



lN/2] 



[(iV-l)/2] 



X21+1 



1=0 



1=0 



where 



-n/A 



Xn 



X 



X 



g jJ^sinh7(gj; 
and 5 is the Kronecker symbol. 



Fl[q] coth(-^7fe)+A"M5(E^^' 27r 



i=i 



5 Scaling limit 

An important stage in the study of the two-dimensional Ising model is the analysis of 
its scaling limit [316]. A new effect with respect to the case of the infinite plane is that 
the "cylindrical parameters" Vi which tend to zero for fixed distance from the 

critical point (A^ — >■ cxo, sinh2A' — 1 = const 7^ 0), do not vanish in the scaling limit 
on the cylinder. Although the formulas for correlation functions obviously become more 
involved, it is possible to generalize several results obtained for the infinite plane to the 
case of the isotropic Ising model on the cylinder [7]. For instance, it was shown [8] that 
the Ising model correlation functions on the cylinder satisfy some integrable equations 
generalizing Painleve III and Painleve V equations obtained in the planar case. When 
considering the scaling limit of the anisotropic model on the cylinder, we do not encounter 
much difficulty provided we take into account some subtleties in the definition of scaling 
variables. 

Since in this case there are two parameters K^, Ky, instead of a critical point we have 
the critical line sinh 2Kx sinh 2Ky — 1 = 0. One has certain freedom in the definition of 
the distance from the critical line, which can be used to ensure that the results obtained 
in the anisotropic case in the scaling limit coincide with those obtained for the isotropic 
{Kx = Ky) model. It turns out that this condition is satisfied with the following choice 
of scaling variables: 

Ir^l — >• 00, \ry\ — i> 00, N ^ 00, 7(0) 0, 

7(0)rj. = X = const, 7 (0)ry sinh 27^'^. = ?/ = const, 7(0)A^ sinh 2/^3, = /3 = const. (40) 

Because of the presence of the exponential factors e"'^'"''''*^'^-'"^*''""'^ and e"^'^^'^-* in the corre- 
sponding sums and integrals, we can suppose that g <^ 1 and, therefore. 
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In order to find the scaling limit asymptotics of the cylindrical parameters A, C,t and r]{q), 
it is convenient to use integral representations 



A ^ = — dp Incoth 



ri{q) 



6 



1 

vr 



Ar2 



dp 



cosp — e 



-7(g) 



cosh 7(g) — cosp 
dpdq^'ip) 7'(g) 



In coth 



Nl{p) 



J sinhA^7(p) sinhA^7(g) 


equivalent to the formulas (fTSi) . (120|) . (l25l) . Denoting 



In 



sin(j9 + q)/2 



sin(p — q)/2 



uj{q) = ^/l + q^ 

and using the above relations, we obtain the answer for the scaled correlation function of 
the 2D anisotropic Ising model on the cylinder: 



n=0 

00 

(a(0,0)a(r.,r,))(+) = ^tI/^) e-N/^('3) ^ 1/, /5) 



n=0 



gn{x,y,p) 



nip''' 



En 



-|^l'^('?j)+«?/<?j-»?{9j:/3) 



Here, the summation is performed over integer / 



uj{qj 



(41) 
(42) 

(43) 



[Z] /i=— 00 l„=-co 



the quasimomenta take on the bosonic values qj 
and form factors are given by 



2itL 



, and the cylindrical parameters 



dp Incoth 







- I dp .'"^^^-^ . In coth .^^^^^ 







P"^ f dp dqui'ip) uj'{q) 

2tx'^ J sinh/3ci;(p) smh.puj{q) 


13 



In 



p + q 



p-q 



l<i<j<n ^^1' 

We see that the spin-spin correlation function on the cyhnder in the scahng hmit defined 
by fHOj) is given by the same functions of renormahzed coordinates as in the isotropic 
model. Therefore, all results of [8] (determinant representations of the correlation fuctions, 
differential equations) apply to the anisotropic well. 



6 Spin matrix elements 

In order to calculate multipoint correlation functions on the cylinder and torus, it is neces- 
sary (and sufficient) to have formulas not only for the form factors \ns{^(^{.^i^)\'Pi ■ ■ - Vl) 
but also for all other spin matrix elements. For instance, the two-point correlation func- 
tion in the Ising model on the periodic lattice of size M x in the ferromagnetic region 
of parameters can be written as follows: 



(a(0,0)fr(r,,r^)) 



Tr {(t(0, 0)T'^-a(0, ry)T^'^~''- ] 



Tr r 



M 



2K 2L 

[N/2+1] [N/2+1] T. Qj-iry E Pj 

2^ 2^ 2^ 2^ \Ns\qi- ■ ■q2K\(y{^,^)\pi- ■ ■P2l)r\ 

K=Q L=0 qi-..q2K P1---P2L 



{2K)\{2L)\ 



2K / 2L N 2K / 2L N 

r -(M-r^) E l{lj)-r. (a-i+ E i{Pj)) -r. E l{q3)-{M-r^) U-^+ E 7(p,)) 1 / f f ^ ^ 

xje +e I / Tr f — 

where Aq denotes the largest transfer matrix eigenvalue, which corresponds to the eigen- 
vector |0)Ars- 

In [9j, we have proposed a formula for arbitrary spin matrix elements of the isotropic 
Ising model on a finite periodic lattice. This formula can be easily generalized to the 
anisotropic case, using two hints from the above: 

• assume that all the difference due to anisotropy consists in the redefinition of the 
function 7(g) in the formula for spin matrix element, and in the appearance of the 

ty{\ — tj. 
tJl - t? 



factor -j-j- 1- to some power (see, e. g., the formula fl36l) ^ 



• also assume that in the scaling limit pol) all multipoint correlation functions are 
given by the same formulas |71 1^ as in the isotropic case. 

Omitting the calculations, we present only the final expression for the spin matrix element, 
which follows from the above assumptions: 

\Ns{qi ■ ■ ■ qK\cr{^, 0)|pi . . .Pl)r\ = ^^t\\ .r . . f \ \ AT U ( \ ^ 
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tji-ti 



K 

n 

Kj 



sm 



2 qj-qj 



sinh 



2 7fe)+7(gj) 
2 



n 



sm 



2 Pi-Pj 



„;„u2 7(P»)+7(Pj) 

bimi 2 i<j<i<' 

i<i<i 



n 



sinh 



2 7(gi)+7(Pj) 



sm 



2 g—Pj 



(44) 



As in the isotropic model case, this representation was verified for finite-row Ising chains 
with N = 1,2,3,4. The scaling limit of the expression fH4|) obviously coincides with the 
classical result [11] in the limit of continuous infinite plane, and with the corresponding 
results on the continuous infinite cylinder [TO] . 



7 Discussion 



Lattice systems with cylindrical geometry recently became important for nanoelectronics. 
However, the elementary cell in these systems is hexagonal. The star-triangle transfor- 
mation brings the Ising model defined on such a lattice to the Ising model on a triangular 
lattice. The simple form of the expressions obtained in the present paper gives a hope 
that similar results can also be obtained in the latter case. 

Exact expressions for the correlation functions are also known for several ID quantum 
systems, in particular, for the XXZ Heisenberg spin chain [12]. In these models, matrix 
elements of local operators on the hamiltonian eigenstates were calculated using the Bethe 
ansatz. One can try to use our results to analyze the behaviour of such integrable one- 
dimensional models in a finite volume and/or at non-zero temperature. 

We thank S. Z. Pakuliak, V. N. Roubtsov and V. N. Shadura for their interest to 
this work and useful discussions. This work was supported by the INTAS program under 
grant No. 03513350. 
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